Project supported by the Key Science and Technology Research Project of Henan Province, China (Grant Nos. 162102210164 and 1721023100107) and the Natural Science Foundation of Henan Educational Committee, China (Grant No. 17A140002).
Abstract
A dualband terahertz (THz) absorber including periodically distributed cross-shaped graphene arrays and a gold layer spaced by a thin dielectric layer is investigated. Numerical results reveal that the THz absorber displays two perfect absorption peaks. To elucidate the resonant behavior, the LC model is introduced to fit the spectra. Moreover, the strength and linewidth of the absorption peak can be effectively tuned with structural parameters and the relaxation time of graphene. Owing to its rotational symmetry, this THz absorber exhibits polarization insensitivity. Our designed absorber is a promising candidate in applications of tunable optical sensors and optical filters.
The metamaterial perfect absorber (MPA) was experimentally discovered by Landy et al. in 2008 and has attracted extensive research interest.[1] The MPA is a kind of artificial plasmonic nanostructure that is usually composed of a periodic metallic structure on a ground plane and spaced by a dielectric layer. By carefully and reasonably designing its physical size and material parameters, it can couple with the incident electromagnetic wave to achieve perfect absorption of the electromagnetic wave in the GHz frequency band. Since their discovery, MPAs have experienced rapid development due to their applications in thermal, biomedical, chemical sensing, and other fields. The terahertz (THz) wave is located in a special area of transition from electronics to photonics. Through scale optimization, new devices based on metamaterials can also produce an effective response in the THz band. Ye et al. numerically demonstrated that a composite structure of a cut-wire array can be used as an effective absorber with a resonant absorption up to 99.9%.[2] Wang et al. designed a mechanically tunable absorber in the THz region. The resonant shift was achieved due to the geometric size of the absorber.[3] Luo presented a thermally dependent multiband THz absorber, which was comprised of a periodic array of closed metallic square ring resonators and metal bars.[4] In general, there are two ways to achieve multiband or broadband absorption of the incident light.[5,6] One method is to design multiple resonators with different sizes and then combine them together.[7,8] The other method is to stack resonators of different sizes on multiple dielectric layers.[9,10] However, the aforementioned structures are difficult to initiate in experiments because of the complicated metallic structure and high cost.[11,12]
Recently, graphene has been introduced into the design of plasmon-induced absorption devices due to its unique optical and electrical properties.[13,14] It exhibits the characteristics of electric field tunability and low propagation loss. In particular, the Fermi energy Ef of graphene can be dynamically tuned by using the external electrostatic gating.[14,15] Rasoul et al. achieved a frequency-tunable MPA based on graphene micro-ribbons on a thick dielectric layer deposited on top of a reflecting metal substrate.[16] Luo et al. showed a periodic array of graphene nanodisks exhibiting 100% light absorption.[17] In addition, the combination of graphene wire and gold cut wire have been studied.[18] However, such a design usually embodies only singleband or narrowband absorption. Khavasi proposed a method to achieve broadband absorption with a periodic array of graphene ribbon, the normalized bandwidth absorption of which can reach up to 1.0.[19] Nonetheless, such proposed structures are polarization-dependent, which hampers their potential applications.
In this work, a cross-shaped graphene THz absorber is proposed. It can achieve THz dualband absorption. The designed structure consists of a periodically distributed cross-shaped graphene layer and a gold layer spaced by the SiO2 dielectric layer. Interestingly, the resonant absorption peaks display polarization independence, namely, the resonant peak position remains the same under different polarization angles. The physical mechanism of the two resonant modes is revealed with the LC circuit model. The low-frequency resonant peak f1 results from a typical dipole, and the high-frequency resonant peak f2 shows a pair of opposite phase dipoles. Benefitting from the unique properties of graphene, the absorption of the proposed structure can be dynamically tuned through the Fermi energy and the intrinsic relaxation time via electrostatic gating instead of refabricating the structures. Moreover, the figure of merit (FOM) of the absorber is 15.35, which is higher than that of the same type of previous absorbers (14.55 in Ref. [20]). Due to its special tunable characteristics, our designed graphene structure is a promising candidate for applications in tunable optical filters and sensors.
2. Structure design and simulation method
The designed structure of the periodic cross-shaped graphene THz absorber is presented. Figure 1(a) displays the unit cell of graphene nanostructures on the SiO2 dielectric spacer. Its dielectric coefficient is 3.9 and the thickness is . a denotes the width of the cross-shaped graphene, which is , and b denotes the length of the cross-shaped graphene, which is . It is periodic in the x and y directions with periodicity , as shown in Fig. 1(b). A linearly polarized plane wave light is normally incident along the z-axis and the angle of the polarization direction with respect to the x-axis is defined as θ. The ground plate is gold with a conductivity , which is a perfect reflection in the frequency domain of interest. Both the transmission spectra and the field distribution are calculated by using COMSOL Multiphysics. It is a commercial software based on the finite element method. It subdivides a large system into smaller, simpler parts that are called finite elements. The simple equations that model these finite elements are then assembled into a larger system of equations that models the entire problem.[21] The solution of Maxwell’s equations is used to calculate the electric and magnetic field components of each point in the grid space. With the continuity of the tangential electric field at each boundary, the parameters of the simulated space electromagnetism are all given in the space grid. In our calculation, the computation domain consists of only a unit cell with the perfectly matched layer applied to the boundaries along the propagation direction z and periodic conditions employed to the four lateral boundaries in the x and y directions.
Fig. 1. Schematic view of the graphene-based dualband THz absorber. The geometrical parameters of the proposed structure are as follows: , , , the thickness of SiO2 is , and the thickness of gold is .
To carry out the simulation, graphene is modeled as a conductive surface.[18,22,23] The transition boundary condition is used for graphene, and its thickness is set as 1 nm. The conductivity of graphene can be derived within the random-phase approximation in the local limit:[17,24,25]
Here, kB is the Boltzmann constant, T is the temperature, ω is the frequency of light, τ is the carrier relaxation lifetime, and Ef is the Fermi energy. In the calculation, our method mainly focuses on the material effect of the two-dimensional flat surface while ignoring that in the out-of-plane direction. At room temperature and low THz frequency (, , the in-plane conductivity of the graphene can be represented by a Drude model:[26,27]
The intrinsic relaxation time is expressed as , where μ is the measured carrier mobility, and m/s is the Fermi velocity.
3. Results and discussions
The calculated absorption spectra of the cross-shaped graphene-based structure is illustrated in Fig. 2. Two perfect absorption peaks appear, corresponding to the absorption value of 96% at 6.25 THz and 97% at 14.25 THz with Ef = 1.0 eV and τ = 1.0 ps. To fully determine the physical origin of the dualband THz absorber, the normalized electric field distribution and the z-component electric field distribution at the absorption peak are illustrated in Fig. 3. Figure 3(a) demonstrates that the electric field concentrates mostly on both ends of the vertical graphene ribbon, corresponding to the absorption spectra of the mode f1 at 6.25 THz. Moreover, we can find the electric dipole resonance from the z-component electric field distribution in Fig. 3(c). Figure 3(b) shows that the electric field is mostly located around the four corners of the cross-shaped graphene, and Figure 3(d) shows a pair of opposite phase dipoles that are similar to a quadrupole. Specifically, the electric fields are mainly distributed on the vertex of the graphene ribbon in the horizontal direction.
Fig. 3. The normalized electric field distribution of the dualband THz absorber at the absorption peaks (a) 6.25 THz and (b) 14.25 THz, and the z-component electric field distribution at the absorption peaks (c) 6.25 THz and (d) 14.25 THz.
In order to get a clear distribution of the resonant absorption peak position, the calculated absorption spectra of the proposed structure under different aspect ratios a:b (a is taken as a constant) are shown in Fig. 4. With the decrease in length b, the peak position of the mode f1 has a significant blue shift in Fig. 4(a). The resonant frequency of mode f2 also increases with the decrease in length b but the degree of blue shift is much smaller than that of mode f1 in Fig. 4(b). Employing the LC circuit model, we can gain further insight into the resonant peak position change of the f1 mode. The resonant frequency can be obtained with the following equation:
where C is the capacitance and it can be approximately described by a capacitor formula .[28]L is the inductance of the structure, which can be given as . From Eq. (3), it is clearly found that a smaller b leads to an increase in the resonant frequency of mode f1. This agrees well with our simulation results in Fig. 4(b), namely, the resonant absorption peak shows a significant blue shift with an increase in the aspect ratio . Because the change of the resonant mode f2 mainly depends on the length of a, the peak position of mode f2 is not very different from that of mode f1.
Fig. 4. (a) Absorption spectra of the dualband THz absorber with different aspect ratios a:b=0.3, 0.4, and 0.5. (b) The change in peak position as a function of the aspect ratio.
The conductivity of the graphene layer depends to a large extent on the Fermi energy Ef, which can be tuned through the external electrostatic voltage or chemical doping. The absorption spectra of our proposed structure at different values of Ef are shown in Fig. 5. It is found that the resonant peaks move to the high-frequency region with the increase in Fermi energy. As we know, the resonant wavelength of the graphene can be given as , where α and β are constants related to the geometry parameters and surrounding dielectric properties, respectively, and ng is the effective index of the graphene.[23,29] The effective index of the graphene layer decreases when the Fermi level changes from 0.8 eV to 1.0 eV. This results in a decrease in the resonant wavelength. Thus, an increase in the resonant frequency occurs with a larger Ef in Fig. 5.
Fig. 5. Absorption spectra of our proposed dualband THz absorber with different values of Fermi level EF.
Relaxation time τ is usually used to describe the gradual return from a certain state to an equilibrium state. It can be expressed as ,[30] where μ is the measured carrier mobility, m/s is the Fermi velocity, and Ef denotes the Fermi energy. Therefore, τ can be dynamically tuned by changing the Fermi energy or mobility of graphene.[31] The Fermi energy of graphene can be tuned via electrical gating or chemical doping. For gate tuning of the chemical potential, the ion gel gating method can be adopted using the field-effect transistor structure with the graphene supported by fused silica.[32–34] Thus, the electrical tunability of relaxation time in graphene can be realized in the following procedure. Firstly, a uniform graphene sheet obtained with the chemical vapor deposition method is transferred onto an ITO-coated silica substrate and then patterned into cross-shaped arrays by using electron beam lithography and oxygen plasma etching.[35,36] After patterning an array of graphene, we cover it with ion gel. A bias potential is applied between the ITO and a gold contact to electrically dope the graphene through the ion gel, which allows us to controllably tune the chemical potential over a sufficiently large range.
The absorption spectra of our proposed structure at τ in graphene elements with Ef = 1.0 eV are shown in Fig. 6(a). The varied intrinsic relaxation time in graphene can affect the absorption strength and linewidth of the absorption peak. The carrier mobility will be significantly enhanced when organic molecules are placed on graphene,[37] leading to an increase in relaxation time. The changes in the absorption and full width at half maximum (FWHM) at the resonant frequency as a function of the relaxation time τ are shown in Fig. 6(b). When the relaxation time ranges from 0.25 ps to 2.0 ps, the FWHM at these two resonance modes gradually decreases. In order to obtain the best absorption spectrum of graphene, a parameter Q is used to describe the relationship between the absorption intensity and the relaxation time. Q is defined as the intensity of absorption A divided by the FWHM, namely, Q=A/FWHM. The detailed performance of FWHM and A with different relaxation times τ is summarized in Table 1. It can be observed that the variation of FWHM and A is not consistent with the increases in relaxation time. The absorption intensity firstly increases and then decreases with increasing relaxation time. The absorbance and FWHM can be comprehensively evaluated by the parameter Q. The larger Q indicates that the absorption intensity and FWHM of the absorption spectrum have ideal values at the same time. When τ =1.0 ps, parameter Q reaches a maximum of 2.67. Therefore, we choose τ =1.0 ps to obtain the absorption spectra with higher absorption intensity and narrower FWHM.
Fig. 6. (a) Absorption spectra of our proposed dualband THz absorber with different relaxation times τ. (b) The FWHM at the resonant frequency as a function of the relaxation time.
Table 1.
Table 1.
Table 1.
Calculated values of parameters FWHM, A, and Q.
.
τ/ps
0.25
0.5
1
2
A
FWHM
Q
A
FWHM
Q
A
FWHM
Q
A
FWHM
Q
f1
0.79
0.99
0.80
0.97
0.62
1.56
0.95
0.43
2.21
0.75
0.35
2.14
f2
0.81
0.91
0.89
0.98
0.55
1.78
0.96
0.36
2.67
0.75
0.32
2.34
Table 1.
Calculated values of parameters FWHM, A, and Q.
.
As an application of sensing, the absorption peak is sensitive to the variation of a nearby dielectric medium. Figure 7(a) illustrates the absorption spectra with different refractive indices of the substrate. When n is larger, the position of the absorption peak shows an obvious red shift. The refractive index sensitivity of the resonant modes can be qualified with the FOM. It can be expressed as follows:[38]
where indicates the frequency shift per refractive index unit, which can be obtained from the fitted slope. Combined with the results obtained from Fig. 7(b), it is found that the FOM of f2 can be up to 15.35, which is higher than that of the current best FOM absorber (14.55 in Ref. [39]). Therefore, the f2 absorption peak with a high FOM can be employed in sensing applications.
Fig. 7. (a) Absorption spectra of our proposed dualband THz absorber with different values of the refractive index n. (b) The relationship between spectral shift of the absorption peak and the refractive index of the substrate.
The transmission spectra of the dualband THz structure under different polarization angles θ are investigated in Fig. 8. Electric field polarization angle θ is defined as the degree off the y axis. Here, θ is taken as 0°, 45°, and 90°, respectively. It can be clearly observed that the resonant absorption peaks remain stable under different polarization angles, embodying polarization independence, which is normally due to the rotational symmetry of the unit structure.
Fig. 8. Absorption spectra of our proposed dualband THz absorber under different incident angles.
4. Conclusions
In summary, we propose and numerically investigate a dualband THz absorber constructed with a graphene layer and gold film spaced by a dielectric layer. Numerical simulations indicate that two perfect absorption peaks at 6.25 THz and 14.25 THz can be realized. A simplified LC model and the electric field distribution are introduced to elucidate the resonant behavior. Our proposed absorber also embodies the polarization independency. Moreover, the absorption peaks can be dynamically tuned through the Fermi energy and the intrinsic relaxation time. With the performed sensitivity measurements, the FOM in the graphene-based structure can reach up to 15.35, which is higher than that of the same type of previous absorbers. Thus, our absorber can be used in many promising sensing applications.